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ABSTRACT 

Recently, the gamma-ray telescopes AGILE and Fermi observed several middle-aged 
supernova remnants (SNRs) interacting with molecular clouds. A plausible emission 
mechanism of the gamma rays is the decay of neutral pions produced by cosmic ray 
(CR) nuclei (hadronic processes). However, observations do not rule out contributions 
from bremsstrahlung emission due to CR electrons. TeV gamma-ray telescopes also 
observed many SNRs and discovered many unidentified sources. It is still unclear 
whether the TeV gamma-ray emission is produced via leptonic processes or hadronic 
processes. In this Letter, we propose that annihilation emission of secondary positrons 
produced by CR nuclei is a diagnostic tool of the hadronic processes. We investigate 
MeV emissions from secondary positrons and electrons produced by CR protons in 
molecular clouds. The annihilation emission of the secondary positrons from SNRs 
can be robustly estimated from the observed gamma-ray flux. The expected flux of 
the annihilation line from SNRs observed by AGILE and Fermi is sufficient for the 
future Advanced Compton Telescope to detect. Moreover, synchrotron emission from 
secondary positrons and electrons and bremsstrahlung emission from CR protons can 
be also observed by the future X-ray telescope NuSTAR and ASTRO-H. 
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1 INTRODUCTION 

The origin of cosmic rays (CRs) is a longstanding problem 
in astrophysics. Supernova remnants (SNRs) are thought 
to be the origin of Galactic CRs. The most popular SNR 
acceleration mechanism i s diffusive shock acceleration (e.g. 
iBlandford fc Eichledll987h . In fact, Fermi and AGILE show 
that middle-aged SNRs (~ 10 4 yrs old) interacting with 
molecular clouds emit GeV gamma rays with broken power 
law spectra (e g.lAbdo et alj|2009l . |2010J ; iTavani et aLlbuTOl ; 
iGiuliani et al.ll201ll) (Hereafter, these references are referred 
to as OBS). Considering a stellar wind before the super- 
nova explosion, the molecular cloud has been swept away 
by the stellar wind. The inner radius of the molecular cloud 
becomes typically about a few tens of parsecs. Then, the 
SNR collides with the molecular cloud at typically 10 4 yrs 
later and the broken power law spectra of GeV gamma 
rays can be interpreted as the inelastic collision between 
molecular clouds and CR nuc lei that have escaped from 
the SN R (hadronic processes) l|Ohira, Murase fc Yamazakil 
l2011ah . However, observations do not rule out contributions 
from bremsstrahlung emission due to CR electrons (lep- 
tonic processes). SNRs have been also observed by TeV 
gamma-ray telescopes, and it is still unclear whether the 
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TeV gamma-ray emission is produced via inverse Compton 
scattering of CR electrons (leptonic processes) or hadronic 
processes. There are many TeV-unID sources in our galaxy, 
and their emission mechanisms are also still unclear. Old 
or midd le-aged SNRs are the candidates of the TeV-unID 
sources l|Yamazaki et al.ll2006l ; lOhira et al.ll2011bh . Thus, it 
is crucial to identify the emission mechanism of gamma rays. 
Positrons and neutrinos are produced in hadronic processes. 
On the other hand, they are not produced in leptonic pro- 
cessesQ Therefore, direct or indirect detections of these par- 
ticles will enable us to identify the emission mechanism in 
gamma-ray sources. 

Recent ly, the CR positron excess has been provided by 
PAMELA (|Adriani et alj|2008h . Although the origin of this 
excess has been invest i gated in many theoret i cal studies 
l|Kashivama et al.l 1201 ll ; Kawanaka et all l20ld ; llokal l20ld . 
and references therein), it still remains unclear. SNRs have 
been discussed as t h e cand i dates of the positron sources 
l|Fuiita et all l2009al : iBlasil 120091 ; IShaviv, Nakar fc Piranl 



Strictly speaking, leptonic processes can produce positrons via 
e~ +7 — > e~ +e — +e + but this reaction rate is suppressed by the 
fine structure constant compared with the gamma-ray production 
rate of inverse Compton scattering. 
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Figure 1. Timcscalcs of interaction for CRs in a molecular cloud 
with the number density n = 300cm -3 , the size R c = 30 pc and 
the magnetic field B = 30 fiG. The blue and green lines show 
the cooling times of secondary positrons (and electrons) and CR 
protons, respectively. The red line shows the annihilation time 
of positrons. The black solid line shows the diffusion time with 
X = 0.01. 



l2009l ; lBiermann et al.ll2009l ; lMertsch fc Sarkarll2009T ). Hence, 
it is important to observe the positron production at SNRs. 

Gamma-ray telescopes showed that the photon flux 
above 100 MeV from SNRs with an age of about 
10V is F >100 m c v = /~ 0MeV d£(dF/d£) = 10 " 7 ' 
10" 6 photons -1 cm -2 (OBS). If the gamma rays have a 
hadronic origin, about the same numbers of positrons are 
produced. When those positrons lose their energy suffi- 
ciently, about 80 percent of them would make positro- 
niums with ambient electrons, and one fourth of them 
woul d decay into two pho tons with the energy of 511 keV 
(e.g. iPrantzos etHI l201Ch . This means that we can ex- 
pect the 511 keV photon flux of the order of 10 -8 — 
10~ 7 photons -1 cm -2 . This flux is sufficiently high for five- 
years observations of the future Advanced Compton Tele- 
scope (ACT) to detect. 

In this Letter, considering the interaction between 
CR protons and molecular clouds, we estimate anni- 
hilation emission of secondary positrons produced by 
the CR protons. We partly refer to analyses of previ- 
ous works concerning the 511 keV lin e from the Galac- 
tic center llAgaronvan fc Atovanl Il98ll ; iBeacom fc Yiiksell 
l200d:ISizun. Casse fc Schannd 2006) and the Galactic plane 
|Steckerl Il967l . Il969l l. and we apply their treatments to 
middle-aged SNRs interacting with molecular clouds. Re- 
cent review o f pos itron annihilation can be found in 
IPrantzos et~ai1 l|201fj ). 

We first provide some timescales for CR protons, sec- 
ondary positrons and electrons inside a molecular cloud 
(Section [2}. We then solve energy spectra of the secondary 
positrons and electrons (Section|3J and calculate the annihi- 
lation line flux as well as the continuum spectrum (Section 
UJ. Section [5] is devoted to the discussion. 



2 RELEVANT TIMESCALES 

In this section, we briefly summarize relevant timescales of 
CRs in a molecular cloud. Physical quantities of molecular 
clouds associated with gamma-ray sources have not been un- 
derstood in detail. Observations suggest that the gas number 
density is about n = 10 2 — 10 3 cm -3 , the size is of the order 
of 10 pc (OBS), and th e magnetic field in mo lecular clouds is 
about B = 1 - 100 /iG l|Crutcher et alj|20ld ). In this Letter, 
for the molecular cloud we adopt values of the gas number 
density n = 300 cm -3 , the size R c = 30 pc and the magnetic 
field B = 30 /iG. We consider ionization and inelastic col- 
lisions with nuclei as cooling processes of CR protons. On 
the other hand, we include ionization, bremsstrahlung, syn- 
chrotron emission, and inverse Compton scattering with the 
cosmic microwave background (CMB) as cooling processes 
of secondary positrons and electrons. For SNRs observed by 
Fermi and AGILE, SNRs interact with molecular clouds and 
the expansion of the SNRs is strongly decelerated. There- 
fore, we here do not consider the adiabatic cooling. 

We here define r coo i = E/E as a cooling time, where 
E and E are the kinetic energy and the energy loss rate, 
re spectively. We adopt an expressio n of E for protons given 
in iMannheim fc Schlickeise"rl (Il994 ), and that f or po sitrons 
and electrons given in IStrong fc Moskalenkol |l998). The 
timescale of direct annihilation of positrons and electrons 
is represented by 



1 



(1) 



n^a a (E)v(E) ' 

where n± , a a and v are the number density of positrons or 
electrons, the cross section of the annihilation I Diradll930r i 
and the velocity of secondary electrons or positrons, respec- 
tively. The electron density, n_, is dominated by thermal 
electrons and the positron density, n+, is very small, so that 
we can neglect the annihilation of secondary electrons. 
The escape time through diffusion is written by 



ME) = 



Br 



W(E) 
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where D(E) is the diffusion coefficient. It is notable that 
so far the diffusion coefficient around an SNR has not been 
understood well. Thus we assume the diffusion coefficient in 
molecular clouds to be 



D{E) 



10" 
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where c is the velocity of light. He re y = 1 might be rea - 
sonable as the Galactic mean value l)Berezinskii et al1ll990h , 
but y = 0.01 should be pr e ferred in surroundi n gs of SNRs 
llFuiita et al l l2009bl . |20ld : iTorres et al.1 l20ld ; iLi fc Chenl 
2010j), so that we use y = 0.01 in this Letter. 

In Fig. Q] we show those relevant timescales. The blue 
line shows the cooling time of secondary electrons and 
positrons where the synchrotron cooing dominates above 
100 GeV, the bremsstrahlung cooling dominates from 1 GeV 
to 100 GeV, and the ionization loss dominates below 1 GeV. 
The green line shows the cooling time of CR protons where 
the pion production cooing dominate above 1 GeV and the 
ionization loss dominates below 1 GeV. CR protons (about 
1 GeV) producing most positrons do not cool and escape 
as long as the SNR age is smaller than the cooling time of 
1 GeV protons. Cooling times of ionization, bremsstrahlung 
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Figure 2. Energy spectra of CR protons, secondary positrons 
and electrons for s = 2.8 and nt = 2.8 X 10 14 cm" 3 s. The green, 
red and blue lines show CR protons, secondary positrons and 
electrons spectra, respectively. 



emission and inelastic collision are inversely proportional to 
the gas density. Therefore, spectral evolutions of all CRs be- 
low 1 GeV depends on only the density. The black line shows 
the escape time due to diffusion with \ — 0.01. The escape 
of secondary electrons and positrons is negligible compared 
with the energy loss. 

The typical energy of positr ons produced by the CR 
protons is about E t ~ 100 MeV (|Murphvl Il987h . Then its 
cooling time is about 3 x 10 4 (n/300 cm" 3 ) -1 yr. Therefore, 
SNRs with an age longer than its cooling time can produce 
low energy positrons. Once positrons have cooled to 100 eV, 
they form positroniums through charge exchange. There are 
four possible spin configurations of the positroniums. One 
of them has the total spin (singlet) and the three others 
have the total spin 1 (triplet). The singlet state produces 
511 keV line photons, and the triplet state produces contin- 
uum photons below 511 keV. Because both the singlet and 
the triplet state are produced at the same rate, one fourth 
of positroniums can decay into two 511 keV line photons. 

It is remarkable that the annihilation time is approxi- 
mately 10 times longer than the cooling time of positrons 
at around 10 — 100 MeV. Then while the positron is being 
cooled, approximately 10-20 percent of positrons directly 
annihilate with electrons. They produce continuum photons 
above 511 keV. 



3 ENERGY SPECTRA OF SECONDARY 
POSITRONS AND ELECTRONS 

In this section, to calculate the photon spectrum from sec- 
ondary positrons and electrons produced by CR protons, we 
calculate energy spectra of positrons and electrons, f± (t, E). 
An evolution of CR spectra is described by 



+ 



/ 



Ta(-E) T d (E) 



= Q(t,E) 



(4) 



where Q(t, E) is the source term of CRs, the third and the 
fourth terms of the left hand side are sink terms due to 
annihilation and diffusion, respectively. As shown in Fig. [T] 
the diffusion escape for secondary electrons and positrons is 



negligible, so that we neglect the diffusion escape. Then the 
solution for secondary positrons and electrons, f±, can be 
expressed by 



f±(t,E) = 



1 



\E(E)\ 



E (t,E) 



x exp 



de Q±(t',e) 
dc' 



e r*,±\E\ 
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where t is time for cooling from Eo to e, and denned by 
rB ° de' . f e de' 



f = 



\E(e')\ 



t 



\E(e') 



(6) 



and Eo(t,E) is the initial energy of positrons or electrons 
before they cool to E and defined by 

= t . (7) 



\E(e) 

We can neglect the annihilation loss for secondary electrons 
because the timescale is much longer than the SNR age. 

Using the code pr ovided by iKamae et all l|2006l ); 
iKarlsson fe Kamad (|200cf ). we calculate secondary positrons 
and electrons source spectra, Q±, and the gamma-ray spec- 
trum produced by decaying unstable hadrons such as pions 
and kaons. In addition, for the high-energy electron source 
Q_, we consider knock-on ele ctrons produced by Coulomb 
collis ions with CR protons |Abraham. Brunstein fe Clind 
1966). The knock-on electrons contribute somewhat to the 
continuum emission as bremsstrahlung emission. 

To obtain the source term of secondary positrons 
and electrons Q±(t,E), we need to calculate the spectral 
evolution of CR protons. CR protons are injected from 
SNRs to molecular clouds in an energy-dependent way 
(jOhira. Murase fe Yamazakil feoiOV Moreover, the diffusion 
escape from molecular clouds should be considered for CR 
protons above 10 GeV (see Fig. [!}. However, we do not 
need to calculate the precise spectrum of CR protons above 
1 GeV because thanks to AGILE and Fermi, we have al- 
ready known spectra of CR protons above 1 GeV. That is, 
we do not need to calculate the spectral evolution due to 
the diffusion escape for CR protons. On the other hand, we 
have to calculate the evolution of the spectrum of CR pro- 
tons below 1 GeV because gamma-ray spectra do not tell us 
the spectrum of CR protons below 1 GeV. Furthermore, the 
injection of CR protons is thou ght to stop when SN R shock 
collides with molecular clouds (|Ohira et al.ll2011ah . that is, 
the injection of CR protons stopped of the order of 10 4 yrs 
ago. Therefore, we assume that CR protons are injected as a 
delta function in time and an effective steep-spectrum, Qcr, 
instead of neglecting the diffusion escape 

Qcn(t,E) = qcn(E)S(t) , (8) 

where qcR(E) is 

qcn(E) oc (E + rn p c 2 ) {fi(£+2m p c 2 )}^ , (9) 

where m p is the proton mass. This expression gives an in- 
jection term of p~", where p is the momentum of the CR 
proton. Observations show that spectra of CR protons have 
broken power-law forms and the spectral index above the 
break is s = 2.7 - 2.9 except for W44 (OBS). In this Letter, 
we adopt the single power law with s — 2.8 (This assumption 
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does not affect the flux of the annihilation line significantly) . 
Then, the solution to equation Q for CR protons, /cr, can 
be expressed by 

fc^t,E)= E{E ^ )) q ME (t,E)) . (10) 
E(E) 

As mentioned in previous section, the cooling of CRs be- 
low 1 GeV depends on only the density. Hence, the solution 
(equation (|10|l 'l depends on only nt. 

To calculate source terms of secondary electrons and 
positrons, Q±(t',E), in equation ([5]), we approximately use 
the present spectrum of CR protons, fcn(t,E), by chang- 
ing Q±(t', E) into Q±(t,E) in equation ([5]). This is because 
CR protons (about 1 GeV) producing most positrons do 
not cool and escape for SNRs observed by Fermi and AG- 
ILE. The solution for the secondary positrons and electrons 
(equation ([5}) depends on only nt below 100 GeV. Here- 
after, we use nt as the parameter to describe the system 
(for example, nt = 2.8 x 10 14 cm -3 s for n = 300 cm -3 and 
t — 3 x 10 4 yr). After the cooling time of CR protons above 
1 GeV (nt > 1.9 x 10 15 cm -3 s), all CRs have already cooled 
and we do not expect any emission. 

In Fig. [2] we show the energy spectra of the CR pro- 
tons, secondary positrons and electrons given in equations 
© and |10]), where s = 2.8, nt = 2.8 x 10 14 cm" 3 s. For 
CR protons (green line) , the spectrum below 100 MeV is 
modified from the initial spectrum qcR.(E) by ionization. 
The electron spectrum (blue lines) is dominated by knock- 
on electrons below 10 MeV. 

A steady-state solution is obtained by changing Eq 
to infinity in equation Jj>). The steady-state spectrum of 
positrons below the typical energy of the source term E t is 
approximately obtained by using the following approxima- 
tion. 



.(E) = q+5(E - E t 



(11) 



where q+ is the total production number of the secondary 
positrons per unit time. Then, the steady-state spectrum 
below E t is given by 



ME) 



\E(E) 



exp 



Ex 



Ac' 
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(12) 



4 RADIATION SPECTRUM 

In this section, we calculate the radiation spectrum from 
secondary positrons and electrons at nt = 2.8 x 10 14 cm -3 s. 
In this case, primary CR electrons are negligible as 
long as the ratio of primary CR electrons to CR pro- 
tons is K cp < 0.01 (see Fig. [2]). We calculate syn- 
chrotron emission, inverse Compton emission with CMB and 
brems strahlung emission from secondary p ositrons and elec- 
trons (|Strong. Moskalenko fc Reimerll2000h. bremsstrahlun g 



emission from CR protons ( Schuster fc Schlickeiserl 120031) , 
and annihilation emission from secondary positrons (e.g. 
ISizun et al.l I2006T ). Energy spectra of secondary positrons 
and electrons, obtained from equation ([5]), are shown in 

Fig. m 

When positrons cool to ~ 10 2 eV, they form positron- 
iums through charge exchange. The production rate of the 
positroniums, Qp s , is obtained by 




10 



icr 10 
Photon energy [ eV ] 



10' 



Figure 3. Photon spectrum from an SNR with -F>iooMeV = 
10~ 6 photons" 1 cm" 2 , nt = 2.8 x 10 14 cm" 3 s and s = 2.8. The 
black dashed line shows the total spectrum. The red, blue and 
green lines show the annihilation spectrum, the bremsstrahlung 
spectrum from secondary positrons and electrons, and the 
bremsstrahlung spectrum from CR protons, respectively. 




nt [ crri s ] 



Figure 4. The line shows F511 kev/^>i00 MeV as a function of 
nt for s = 2.8. 



Q Ps (t) = \E(100 eV)\f+(t, 100 eV) 



(13) 



One fourth of positroniums decay into two 511 keV photons. 
Thus the photon flux of the 511 keV line, fsiikev, is given 
by 



.(*) 



(14) 



87rd 2 ' 

where d is the source distance. The line width of the anni- 
hilation lin e of positroniums emitted from m olecular clouds 
is 6.4 keV l|Guessoum. Jean fc Gillardll2005l 'l. We here use a 
Gaussian profile with the 6.4 keV width as the line structure. 

Fig. [3] shows the photon spectrum normalized so as 
to make F>ioom c v = 10~ 6 photon s _1 cm -2 , where nt = 
2.8 x 10 14 cm' 3 s and s = 2.8. We again note that this con- 
dition is typical for middle-aged SNRs observed by gamma- 
ray telescopes. Annihilation emission (red line) dominates 
at around the 511 keV range. It is notable that synchrotron 
emission and inverse Compton emission with CMB induced 
by secondary positrons and electrons do not contribute in 
the energy range of Fig. [3] Moreover, synchrotron emission 
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by secondary positrons and electrons depends on the mag- 
netic field, the maximum energy of CR protons and the 
spectral index of CR protons (s). Bremsstrahlung emission 
of CR protons (green line) also depends on s. Note that 
bremsstrahlung emission of CR proto ns can be observed by 
the future X-ra y telescope, NuSTAR ilHarrison et al J 1200a ) 
and ASTRO-H (|Takahashi et alj|2010h . 

For the estimation of the annihilation line, we can ne- 
glect the escape loss as long as r coo i < ra for E < 100 MeV. 
Moreover, most positrons are cooled by the ionization loss as 
long as B < 1 mG (n/300 cm 3 ) 1 ' 2 , so that the magnetic field 
(B), the size of molecular clouds (R c ) and the normalization 
of the diffusion coefficient (x) are not important. Hence, 
the ambiguity of the annihilation line is only nt. Fig. [4] 
shows the flux ratio, F 511 k c v/-F>ioo MeV, as a function of nt, 
where s = 2.8. The ratio does not depend on the number 
of CRs and the distance. The flux ratio, F511 kev/-F>ioo MeV, 
becomes constant after nt = 2.8 x 10 14 cm -3 s because al- 
most all positrons have cooled to 100 eV by that time. From 
Fig. [4] we expect fsnkcv ~ 10 -7 photon cm -2 s _1 from 
SNRs with nt = 10 14 - 10 15 cm -3 s and with F>iooMeV ~ 
10 6 photon cm -2 s _1 . The expected flux of the annihilation 
line, -FsiikeV, c an be sufficient for five-years observations of 
ACT to detect (|Boggs et alj|2006h . Therefore we can expect 
to detect the annihilation line of secondary positrons pro- 
duced by the CR protons in SNRs observed by AGILE and 
Fermi. 



5 DISCUSSION AND SUMMARY 

In this Letter, we have investigated the MeV emission spec- 
trum due to CR protons from SNRs interacting with molec- 
ular clouds. We found that for typical middle-aged SNRs ob- 
served by AGILE and Fermi, secondary positrons can cool 
to an energy sufficient to make positroniums. We calculated 
annihilation emission of the positrons and other emissions. 
The expected flux of the annihilation line is sufficient for 
the future gamma-ray telescope, such as ACT, to detect. 
Therefore, we propose that annihilation emission from sec- 
ondary positrons is a important tool as a CR probe. More- 
over, synchrotron emission from secondary positrons and 
electrons, and bremsstrahlung emission from CR protons 
can b e also observed by t he future X-ray telescope, NuS- 
TAR (|Harrison et al.ll2005h and ASTRO-H (|Takahashi et al.l 
l2010h . 

All particles with energies smaller than 1 GeV lose their 
energy due to ionization, as shown in Fig. [T] Not only the 
511 ke V line but also the ionization rate is also a probe of CR 
nuclei |Goto et al.ll2008l ; llndriolo et al.ll201(ih . iBecker et all 
(|20 1 if ) proposed that H^~ and Hg~ lime emissions should 
be observed from molecular clouds if there are many CRs. 
Low energy CR protons and kno ck-on electrons might be 
measured by u sing atomic lines (|Gabriel fc Philipg Il979l : 
lTatischefJl2003h . so that atomic lines also become a probe 
of CR nuclei. Quantitative estimations will be addressed in 
future work. Moreover, CR nuclei can produce nuclear ex- 
citation lines by i nelastic collisions and productions of un- 
stable nuclei Ce.g. iNath fc Biermannlll994l ; iTatischefj I2OO3I ; 
ISumma. Elsasser fc Mannheim! 1201 ll ). Therefore, CR com- 
positions around SNRs can be also investigated by the nu- 
clear excitation lines. Recent observations and theoretical 



studies for CR compositio ns are also remarkable (|Ahn et al.l 
|2010| ; lOhira fc IokaHioTll) 

The spectra of secondary positrons and electrons pro- 
duced by CR protons are different from a single power law 
because of their cooling and injection spectra (see Fig. [2}. 
The spectra of secondary positrons and electrons below 
100 MeV would become harder than the CR proton spec- 
trum around 1 GeV. The number of secondary positrons and 
electrons can be larger than that of primary CR electrons 
when the SNR age is larger than about 10 percent of the 
cooling time of CR protons due to the inelastic collision. 
Even when the number of the primary CR electrons is larger 
than that of secondary positrons and electrons, the spectrum 
is affected by the cooling. These cooling effects of positrons 
and electrons might be a reason why radio spectra are dif- 
ferent f rom expected from ga mma-ray spectra observed by 
Fermi |Uchivama et alj|20ich . However, the radio emission 
from SNRs may be originated from inside the SNR, but 
observed gamma rays may be originated from outside the 
SNR. That is, a different component may produce the radio 
emission. We may have to build more detailed models to 
compare the theory and observed date in future. It is an in- 
teresting future work to compare between observations and 
the theory for individual SNRs. 
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